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Abstract 

Effective field theory is applied to finite-density systems with two-particle 
interactions exhibiting an unnaturally large scattering length, such as neutron 
matter. A new organizational scheme, identified for a large number of space- 
time dimensions, allows for convergent analytic calculations in many-body 
systems and is similar to the hole-line expansion of traditional nuclear physics. 
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Effective field theory (EFT) provides a way to systematically improve low-energy observ- 
ables and make predictions for related processes. It has recently been applied to systems 
with an unnaturally large scattering length by adopting an organizational scheme that forces 
an infinite set of diagrams to be summed at leading order [|TJ. 

Extending the systematic language of an EFT to many-body systems would be bene- 
ficial for describing various systems of interest including nuclei, neutron stars, and atomic 
Fermi-Dirac condensates. Although sophisticated approximation techniques already exist 
for describing finite-density systems [@||, error estimates are difficult to quantify and have 
been considered the "holy grail" of many-body physics by Brueckner. 

At finite density, EFTs have given an accurate understanding of dilute systems Q and 
of Fermi surface effects || such as the superconductor pairing gap Also, the Skyrme 

model and EFTs constructed at saturation density || can parameterize finite-density effects 
with only a few parameters. However, connecting bulk properties at finite density, such as 
the equation of state of nuclear matter (an idealized infinite system of protons and neutrons 
interacting only through the strong force), to known free-space interactions between the 
constituent particles requires a consistent organizational scheme (or power counting) at finite 
density. This will also allow many other interesting questions to be investigated, including 
whether nuclear matter binds in the chiral limit 0] and why nuclear matter saturates. 

An EFT treatment of many-body systems automatically dictates the relevance of higher- 
particle contact interactions. For example, three-particle scattering composed of only two- 
body interactions has divergent loop corrections and the addition of an actual three-body 
contact interaction is required to regularize the result [|K|]. For nucleons in the triton, such 



a contact interaction has been shown to be as important as two-body interactions [pT 
Therefore, it is natural to first focus on an EFT description of finite density systems with 
two states, such as neutron matter, since then Pauli's exclusion principle requires three- and 
higher-body contact interactions to be multiplied by powers of momentum and the EFT 
implies these contributions are suppressed. 

For neutrons, restricting the momenta to be below m n and neglecting spin dependence 
will clarify the discussion since pions and other extended structure can be integrated out 
leaving only simple point-like interactions. In addition, the two proposed free-space EFT 
power countings, known as A-counting and Q-counting [Q, yield equivalent results in this 
regime. Therefore, 0(kF/m n ) corrections will be implicit in the following. To make the 
relevant scales transparent, dimensional regularization with power divergence subtraction 
(PDS) [[EJ will be used below to regulate infinite integrals. 

Non-relativistic particles with (spin-averaged) contact interactions are governed by the 
most general lagrangian 

£ eff = At (id t + ^1\n- Co(AtA) 2 + ic 2 [(MVA) 2 + ((VN^)Nf 

+ Cf (A f VA) ■ ((VA^A) + . . . , (1) 



with particle mass M, S- wave interactions Co and C*2, P-wave interaction Cf, and higher 
derivatives and partial waves represented by the dots. Calculating observables using this 
non-renormalizable lagrangian leads to divergences which need to be regulated. One-loop 
scattering (Fig. |l|a) with external relative momentum k = — k 2 ) in the PDS subtraction 
scheme is 
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d 3 q M M, 

-—fr + tk). (2) 



(2tt) 3 k 2 - q 2 + ie 4vr 

This introduces the renormalization scale \i into amplitudes, and so the low-energy constants 
in Eq. (|1|) must also depend on \i to ensure observables are /i-independent. 

Note that p should be of the same order as the characteristic momentum of the process, 
k, so that the relative size of the loop Eq. (fj) is not altered. Both scales should be small 
compared to the scale of underlying physics (m n for neutrons) to ensure the EFT description 
does not break down. Furthermore, two-particle interactions exhibiting a large scattering 
length introduce a small momentum scale into the system (the neutron scattering length 
a s gives |l/a«| ~ 10 MeV). All this can be summarized by defining a generic small scale 
Q = (k,p, l/a s ), and all observables are well approximated by an expansion in Q. 

Matching the on-shell T-matrix to phase-shift data gives expressions for the low-energy 
constants of Eq. ([]]) in terms of scattering observables jn| such as the scattering length a s 
and effective range r e , 

47T/M 1 4vr/M r e 

G ° = H 7T ' C2k = 7 -i / \2 ^7 k ^ • 

p - l/a s Q (p- l/a s y 2 

Since the particle-particle loop Eq. (0) scales like Q, all insertions of Co contribute at leading 
order 1/Q to the T-matrix, whereas insertions of Ci and other contact interactions are higher 
order. This justifies the Lippmann-Schwinger equation for Cq only, which forms a summable 
geometric series, giving |L2] 

T(k) ~ °° ~ 4?r/M (3) 

M«J x _M^^ + lk) i/ aa + ik - W 

Note that the result is /i-independent and produces the first term in the effective range 
expansion by construction. 

A many-body system of fermions with spin-isospin degeneracy g = (2S + 1)(2/+ 1) forms 
a Fermi sea with Fermi momentum kp and density p = gk F /6n 3 . The kinetic energy per 
particle of a free fermion gas is given by the familiar expression 

E kin _ g f d?k k 2 _ 3 k 2 F 



A pJ (2tt) 3 2M 5 2M ' 

Quantum fluctuations allow particle-hole pairs to form, with interactions still given by the 
lagrangian Eq. ([I]). These contributions to the ground-state energy can be represented by 
closed, connected Feynman diagrams and calculated using finite-density Feynman rules Hf§] . 

Two basic building blocks of finite-density diagrams are depicted in Fig. |l|a and b. Since 
the neutron propagator accounts for both particle and hole propagation 

r<h\ — a ( ■ 0{k F -k) \ 

G (kU - Sop [ ko _ k2/2M + ie + ko -k 2 /2M- l e ) ' (5) 

contour integration of the loop momentum shows Fig. |l|a and b actually represent three 
processes: two-particle or two-hole scattering in Fig. |l|a and one-particle, one-hole scattering 
in Fig. [l|b. Pauli blocking is enforced by the following theta functions, which depend on 
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FIG. 1. Scattering with propagators in the same direction (a) and opposite directions (b). 
Also shown are the two simplest contributions to the ground-state energy (c) and (d). 



the center of mass momentum P = + k 2 and arbitrary relative momentum q (with 
q li2 = iP±q), 



%1 - k F ) 6(q 2 - k F ) 
- q - 0(Qi ~ k F ) 0{k F - q 2 ) 
Q- = 9(k F - q x ) 6(k F - q 2 ) 



for particle-particle, particle-hole, and hole-hole scattering respectively. These three pro- 
cesses can be integrated analytically: 
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d 3 q 
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(2tt) 3 P (k - q) 

d 3 q e- 



ie (2tt) 2 
k F 



/(«-s, s) 



(2tt) 3 q 2 - k 2 + ie (2tt) 



f(K,-s) 6(1 -s) 



(6) 
(7) 
(8) 



having defined the dimensionless momenta s = ~P/k F and k, = \a/k F . Particle-particle 
scattering has contributions from momenta up to infinity and so can be written as a reg- 
ularized integration over all momenta, Eq. (0), which brings the /i-dependence in Eq. (|J), 
minus a term constraining at least one propagator to be inside the Fermi surface, which 
replaces the remnant of asymptotic scattering, ik, by the nontrivial function f(n,s). This 
same functional form appears in hole-hole scattering, but in this case both ki and h 2 must 
be inside the Fermi surface, and so P < 2k F (i.e. s < 1). In fact, the functions in Eqs. (||-|D 
have different analytical forms depending on whether s > 1 or s < 1, and the combinations 
which enter ground-state energy calculations are 
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(9) 
(10) 
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Interpreting one contribution to the ground-state energy, Fig. |l]d, in three different ways 
gives the identities: 
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FIG. 2. Four representative contributions to the ground-state energy. 



d 3 kd 3 P 717/ \ fd 3 kd 3 P ~ 



d 3 kd 3 P 




(a0 . ^-^"d-^-^i^it^-Sj ■ < 12 > 

which all integrate to the same result, serving as an excellent check of Eqs. (PHTT|). 

Extending Cj-counting to finite density shows the small parameter is Q = (/cp, /i, l/a s ), so 



all powers of /ci?a s must be summed at leading order . Also, this choice \x ~ &f ensures the 
regularization does not destroy the relative importance of Eq. (|j), since /(«, s)/7r ~ C(l). 
With this scaling, the kinetic energy per particle Eq. is 0(CJ 2 ), and the contribution 
with one C insertion, Fig. p]c, is 

E\ g( g - 1) rfkfP (g_ z lV4 r n2 

The free-space counting Co ~ l/Q makes Eq. ( |l3l) of the same order as the kinetic energy 
and hence nontrivial properties such as nuclear binding can possibly occur at leading order. 
In fact, each additional Co vertex added to Fig. |l]c requires two more propagators and one 
more loop integral, bringing an extra factor of orderQ 




Co £ 3 ~Q° ("J 



and so a// possible insertions of Co, including non-planar graphs, are leading order, making 
calculations intractable. 

However, evaluation of some representative diagrams will show there is still an additional 
organizational scheme which can be exploited. Take, for example, the Feynman diagrams 
in Fig. 0a-d, which are all of the same order in Cj-counting. They are rings of bubbles with 
particle scattering in either the same direction or opposite directions. A loop momentum is 
assigned to each bubble and a propagator Eq. (|5|) to each line. The final loop momentum 
can be chosen to be the center-of-mass momentum P, which traverses the entire ring and 
is shared equally by the legs of each bubble. A contour integration over the propagators 
enforces energy conservation, so not all propagators of a diagram can be particles (or holes). 



The identification ko ~ k /2M is made for power counting purposes |T 
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The contributions to each diagram are: two particle-particle bubbles, one hole-hole bub- 
ble (2pp-lhh) and one particle-particle bubble, two hole- hole bubbles (lpp-2hh) for Fig. ||a; 
3pp-lhh, 2pp-2hh, and lpp-3hh for Fig. 0b; 2ph-lhp and lph-2hp for Fig. 0c; 3ph-lhp, 
2ph-2hp, and lph-3hp for Fig. |2|d. After working through the Feynman rules, simplifying 
the expressions, and dividing by the density p, the energy per particle for each diagram, 
organized according to the aforementioned particle and hole content, is: 



/(«>s) 




d 3 s d 3 n 



(4tt) 2 
96 rd 3 sd 3 K 



6^f(K. S ,s) 2 + e±f(-K-S,s) 



9*f(K-S,sy-3f(K-H,8 



/(-q-s,s) + 9±f(-K-s,s) 



where a factor of 9(1 — s) has been absorbed into /(/c, —s). The factor of "3" in contributions 
to Fig. 0b and d comes from the combinatorics of choosing the second pp or ph pair respec- 
tively and the symmetry factor of Fig. |2|d gives an additional 1/3 suppression compared to 
the other diagrams. 

A simple change of variables shows the iph-jhp and jph-zhp contributions are the same 
for any % and j. Numerically integrating the above expressions and taking p = for clarity 



gives 



A J (2a) 

-) 



2M 

kp 

2M 



(k F a s f (0.0641 + 0.0115) 



(k F a s y (0.0383 + 3 x 0.00247 - 0.000683) 



(g -l)(g-3)^ (k F a s ) 3 (0.0287 + 0.0287) 



(15) 
(16) 
(17) 



k 2 

(g - lf(g - 3)^ (k F a s ) (0.000603 + 3 x 0.000557 + 0.000603) . (18) 



Adding the terms for each of the three-bubble diagrams, Fig. 0a and c, gives results in agree- 
ment with Ref. [|]]. Comparison of individual terms in Eqs. ( p!5| - |T8|) shows each additional 
hole-line is suppressed by about a factor of three. This means every added hole-hole bubble 
is suppressed by almost an order of magnitude. The suppression is the same at nonzero p, 
where k F a s can be taken very large. 

This suppression can be described by a new power counting at finite density which is 
kinematical in nature and will lead to the simplification of many-body calculations. Take 
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the lpp-lhh contribution of Eq. fll2"|) , which using Eq. (^) can be written in terms of one 
relative momentum below the Fermi surface and another above: 



/ 



d 3 Pd 3 kd 3 q 6Z6+ fd 3 Pd 3 kd 3 q e;e + ( k 2 „ /lAlA A , \ 

H k q 'l - — - + o(fc 4 /g ) • (19) 



(2vr) 9 k 2 - q 2 + ie J (2vr) 9 -q 2 + ie\ -q 2 + ie 



The integration region 6^ can be pictured as the overlap of two Fermi spheres (of radius kp) 
with centers separated by the center-of-mass momentum P, whereas the region 6 + is given 
by all momenta excluding these same intersecting spheres [0. In principle, the momenta k 
and q can both be near kp, but the integration measure d 3 P suppresses P ~ 0, which makes 
the approximation g > k appropriate. 

In addition, for P < 2k F , the overlap region of the two spheres is small and suppresses 
additional powers of k 2 in the numerator, since 

£ =-^L + l + 0(l-s), (20) 
4tt -k 2 + ie 2k F K 1 ' v ' 

^^ = ii^ + 0(l- S ) 3 , (21) 

and more generally 6 + /q 2n ~ (9(1 — s)° for every n, and 6^ k 2n ~ (9(1 — s)" +2 . Then, 
working in D dimensions (with the understanding that D = 4 will be taken at the end of 
the calculation) the integral over the center-of-mass momentum gives 

'dss^il-sY-^-j'dss - 2 . (22) 



o D 

for large D. In fact, carrying out the entire calculation in Z)-dimensions shows the suppres- 
sion is actually 1/T> n with T> = 2 D / 2 , which is numerically equivalent to Eq. (^) for D = 4 
but even more convergent in larger number of dimensions. 

From the above, the theta functions which enforce Pauli blocking can be associated with 
powers of IjV: 6 + ~ 0~ ~ 1/V 2 , and Eq. (HD implies 0± ~ 1/D. Therefore, the 

expansion in 1/X> is similar to the traditional hole-line expansion of nuclear physics [0. 
Furthermore, the complicated functions of Eq. (P-pJ]) greatly simplify to 

/(«, s) = 2 + O(I)- 1 ) , /(«, -a) = + 0{V- 3 ) , (23) 

since they always appear with the respective theta functions 6~ and 9 + , and using Eq. (|12|), 

6(a - 1) /(aa, s) = 1 + 0(D- 2 ) . (24) 
3s(s — a) 

Using these simplified forms along with Eqs. (|20| - pTD , the above analytic expressions for the 
diagrams of Fig. |2] can be analytically integrated. The numerical values at each order in 
1/T> are collected in Table |, and the 1/D-expansion quickly converges to the full result in 
most cases, although the contributions to Fig. [|c require further investigation. 

The existence of a power counting in 1/T> allows for systematic calculation of finite- 
density observables such as the ground-state energy. First of all, it is interesting to note 
that at next-to-leading order in the Q-expansion, one insertion of C^k 2 will appear, and the 



7 



TABLE I. Exact factors in contributions to the energy per particle of Eq. (1^-1^) and their 
respective expansion in powers of 1/2?. An overall factor of "3" for the 2pp-2hh and 2ph-2hp terms 
has been left out below for a more direct comparison. The lpp-3hh diagram first appears at order 
1/P 6 with contribution -0.000507. 





Fig. |a 
2pp-lhh lpp-2hh 


Fig. |b 

3pp-lhh 2pp-2hh lpp-3hh 


Fig. |c 
2ph-lhp 


Fig. Id 
3ph-lhp 2ph-2hp 


total 


0.0641 0.0115 


0.0383 0.00247 -0.000683 


0.0287 


0.000603 0.000557 


<D(1/V 2 ) 
0(1/V 3 ) 
0(1/P 4 ) 
0(1/V 5 ) 


0.0860 
-0.0162 

-0.0056 0.0133 
0.0002 -0.0024 


0.0548 
-0.0155 

-0.0032 0.00310 
-0.0007 -0.00164 


0.0283 
0.0057 
-0.0003 


0.000676 0.000676 
0.000371 -0.000230 



integration of k 2 over phase space (8^) brings a factor of 1/T>. Therefore the hierarchy of 
the Q-expansion is already included in the 1/D-expansion. The leading order ground-state 
energy per particle can be written in closed form: 

\ = + + 0(1/V) , (25) 

with the kinetic energy given by Eq. (§) and the interaction energy to this order given by 
all the particle-particle scattering terms (i.e., the npp-lhh bubbles for all n). This is a 
summable geometric series, giving 

Eg = g(g-l) rfPfk 4n/M 
A 2p J (2tt)° k l/a s ->ff(K, S ) ' 1 ) 

which is //-independent, as it should be for any consistent expansion. Actually, to the order 
quoted in l/V, the substitution /(/?, s) — > 2 is valid, which gives the leading-order ground- 
state energy per particle for neutron matter (taking g = 2) 

E 



A 




<D(1/V) . (27) 



For a s > 0, E/A has a pole at kpa s = tt/2 and it would be interesting to explore its 
consequence in atomic systems where the scattering length can be tuned. As for neutrons, 
the scattering length is negative and it is a good approximation to take l/a s — >• 0, to arrive 
at the scale-invariant answer E/A = (4/15)kp/2M. The next order in 1/T> will require a 
single insertion of C2 or Cf and diagrams with three hole-lines. 

A new power counting at finite density has been found which includes the nonperturbative 
free-space power counting of two-nucleon EFTs. The many-body problem simplifies in a 
large number of space-time dimensions and leads to a convergent expansion even for D = 4. 
Including explicit pions and three-body forces in the lagrangian will allow the above analysis 
to be extended to nuclear matter for an exploration of the behavior near saturation density. 
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